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ABSTRACT 


A radiative transfer algorithm in the solar wavelengths 
for the NOAA POES AVHRR is presented for the cloud-free, 
marine atmosphere. This algorithm combines linearized, 
Single-scattering theory with an estimate of bi-directional 
surface reflectance. Phase functions are parameterized 
using an aerosol distribution model and the ratio of 
radiance values measured in channels 1 and 2 of the AVHRR. 
Automated cloud screening and sun glint removal is included. 
Retrieved satellite aerosol optical depth (AOD) is compared 
to surface measured sunphotometer AOD collected during the 
International Global Atmospheric Chemistry (IGAC) Project's 
Second Aerosol Characterization Experiment (ACE-2) from June 
mopetO  JUljy 25, el come The comparison data set has a 
correlation coefficient of 0.88 with a standard error of 
0.02 at both channel 1 and 2 wavelengths. Regional aerosol 
properties are examined with an emphasis on the differences 
between the ACE-1, TARFOX and ACE-2 regions. ACE-1 and ACE- 
2 regions have strong modes at AOD at around 0.1, but ACE-2 
tails toward higher values consistent with urban and dust 
aerosol intrusion. The TARFOX region has a noticeable mode 
at AOD around 0.2, but has significant spread of AOD values 
consistent with the varied aerosol constituents in that 
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I. INTRODUCTION 


Atmospheric aerosols, whether naturally occurring or 
anthropogenic, impact the Earth's energy budget. Incoming 
solar radiation is scattered by these aerosols resulting in 
a net decrease in heating of the Earth's surface. Aerosols 
provide smaller condensation nuclei that increase cloud 
albedo at solar wavelengths, which again reduces heating at 
the Earth's surface. This is exactly the opposite affect of 
the currently highly publicized global warming trend due to 
greenhouse gases. According to Charlson et al. (1992) and 
the Intergovernmental Panel on Climate Change (IPCC) (1996), 
the cooling influence caused by aerosols maybe offsetting 
the greenhouse warming to a substantial degree. Schwartz 
and Andreae (1996) point out that if the aerosol forcing is 
Significant and has negated much of the greenhouse forcing, 
then the resultant increase in global temperatures has come 
from a small residual forcing indicating a greater planetary 
temperature sensitivity. A greater sensitivity may result 
in accelerated global warming in the future. The 
uncertainty of aerosol forcing on the climate needs to be 
reduced to near the levels of uncertainty in greennouse 
forcing in order to compare their weights. One avenue 
towards reducing this uncertainty is through satellite-based 
measurements oof aerosol properties both globally and 
regionally. | 

Characterization of aerosol optical properties is 
equally important for the U.S. Navy. Since aerosols scatter 
and therefore degrade the system signal, this data can be 
used as input for electro-optical/infrared propagation 
models and tactical decision aids. Knowledge of aerosol 


optical properties is also important in the planning and 


design of new weapons and sensors that operate at 
optical/infrared frequencies. 

Ideally, characterization of aerosol optical properties 
would be accomplished through in situ measurements of 
radiative properties and aerosol size distributions. 
Sunphotometers, spectrometers, and radiometers provide very 
accurate pictures of these parameters, but their use are 
limited spatially and temporally. Additionally the 
employment of these tools is expensive and time consuming. 
The use of satellite-based radiometers, as suggested above, 
1S ga vaeable resowroeme that @allows greater  spabiad amd 
temporal coverage. To measure the radiative properties in 
question the National Oceanic and Atmospheric Administration 
(NOAA) Polar Orbiting Environmental Satellite (POES) with 
its Advanced Very High Resolution Radiometer (AVHRR) can 
provide up to two passes per day for local analysis and a 
global picture twice per day. The AVHRR channel 1 (visible) 
and channel 2 (near infrared) possess the appropriate 
spatial and spectral resolution to measure optical radiative 
properties. This two solar channel capability is unique to 
the AVHRR instrument. Work done by Durkee et al (1991), 
Ronault and Durkee (1992) and Brown (1997) explore the use 
of the AVHRR two solar channel capability as a way to 
characterize phase scattering effects of aerosols. Aerosol 
retrieval techniques based on their work are used in this 
study. 

Several experiments over the last few years have 
focused on aerosol properties and have provided ample field 
data for use in closure studies. These experiments include 
the International Global Atmospheric Chemistry (IGAC) 
Project's First Aerosol Characterization Experiment (ACE-1) 
in November 1995, Second Aerosol Characterization Experiment 
(ACE-2) in June 1997, and Tropospheric Aerosol Radiative 
Forcing Observational Experiment (TARFOX) in July 1996. 


These experiments are important because they provide 
information on several different aerosol distributions that 
can represent a majority of the aerosol distributions found 
qlreba lly . ACE-1 focuses on aerosols in the remote marine 
atmosphere; TARFOX focuses on aerosols carried over the 
Western Atlantic Ocean from the United States; ACE-2 focuses 
on anthropogenically-modified aerosols and dust aerosols 
carried over the Eastern Atlantic Ocean from Europe and 
Africa. 


The objectives of this thesis is twofold: 


1) Validate aerosol optical depth retrievals using 


ACE-2 field observations of optical depth, 


2) Conduct a regional analysis based on composites of 
retrieved optical properties during ACE-1, ACE-2, 
and TARFOX. 


Chapter II describes the basic radiative transfer 
theory used in the satellite optical depth retrievals. 
Chapter III outlines the three experiments, the data sets 
and basic instrumentation used. Chapter IV describes the 
aerosol optical depth retrieval procedures and compositing 
procedures. Chapter V presents the results and Chapter VI 


discusses final conclusions and recommendations. 





II. RADIATIVE TRANSFER THEORY 


The mechanisms by which electromagnetic radiation is 
generated and how this radiation interacts with the 
atmosphere is described by radiative transfer theory (see 
Liou, 1980, for a complete description of the theory). The 
remote measurement of atmospheric properties by satellites 
necessitates an understanding of these theories. Salydslis7 
stated, the radiation that is sensed by a= satellite 
radiometer is the sum of the sources: emitted radiation and 
radiation scattered toward the satellite; and the sinks: 
absorbed radiation and the radiation scattered away from the 
satellite. Quantifying these sources and sinks leads to a 
complex equation called the Radiative Transfer Equation 
aT Ee Through assumptions about atmospheric radiative 
properties, the RTE can be simplified to a more usable form. 
In the case of calculating the effects of atmospheric 
aerosols we use a single scatter approximation along with a 
no-emission assumption. In a elllfebrel free, marine 
environment, the short wave radiation sensed by a satellite 
radiometer is due mostly to the scattering of solar 
radiation by the molecular constituents of the atmosphere 
(Rayleigh scattering) and larger suspended aerosols (Mie 
scattering). Short wave here is defined as the visible (0.4 
to 0.7 um) and near infrared (0.7 to 1.0 um) spectrum. 
Corrections for water vapor absorption at these wavelengths 
are on the order of 2-3% and corrections for ozone 
absorption are less than 5%. Although the no-emission 
assumption precludes terrestrial radiation playing a role, 
the ocean may reflect solar radiation. This is considered 
to be small except in certain cases such as sun glint and 
where high concentrations of particulate and phytoplankton 
exist. Contributions from sea surface foam and sub-surface 
reflectance can be accounted EOL with empirical 


measurements; specular reflectance (sun glint) can be 
estimated. The radiative transfer process described here is 
illustrated in figure 2.1. The basic theories used to 
obtain aerosol optical properties in a cloud free, marine 


environment are presented in this chapter. 


Dee OPTICAL CHARACTERISTICS AND RADIATIVE TRANSFER 


The optical characteristics of the atmosphere can be 
inferred through the measurement of their effects on 
radiative processes. Electromagnetic radiation interacts 
with the atmosphere in three basic mechanisms: absorption, 
scattering, and emission. Incoming solar radiation is 
attenuated in proportion to the density and absorbing 
characteristics of the atmosphere. The combination of 
absorption and scattering is termed extinction, and is 


quantified in the extinction coefficient given by: 


TO (m, r)n(r)dx (Zeer) 


ext 


QO 
tI 
on 8 


where r is particle radius, mr* is particle cross-sectional 
area, Q.,,(m,r) is the extinction efficiency factor, m is 
the complex index of refraction, and n(r)is the number of 
particles for a given radius. Changes in the _ size, 
compoSition, and distribution of constituents or suspended 
particles in the atmosphere lead to a change in extinction. 
Therefore, the measurement of extinction can inversely lead 
to the characteristics of the atmosphere's particulate 
dicta butaon . 

The scattering of radiation is the largest contributing 
factor to extinction in the visible and near infrared 
wavelengths. Rayleigh (molecular) scattering can be 


calculated since the atmospheric molecular constituents are 


nearly unaform both spatially and temporally (Durkee et al., 


Poo After the effects of molecular scattering are 
removed from the retrieved radiances, quantification of 
scattering effects due to suspended aerosols are left. The 
effect of scattering due to aerosols iS approximated using — 
Mie theory for spherical particles. For the cloud-free 
marine environment, the extinction coefficient can be 


approximated by the scattering coefficient: 
Osc = | 1’ Qoeas (m0, A,1)(n)dr (oe 
0 


where Q,..,,(m,A4,r) is the scattering efficiency factor, 
representing the ratio of total energy scattered in all 
directions to incident energy. 

Optical depth quantifies the atmospheric extinction 
through vertical integration of the extinction coefficient 


through the atmosphere, as represented by: 


Satellite radiometers measure the radiative properties of 
the entire atmospheric column therefore quantification of 
total column optical depth is the objective of the satellite 
retrieval technique used in this study. 


The above equations show the relationship between 


optical depth and the aerosol particle distribution. The 
extinction coefficient, and therefore the scattering 
Soctmuecnent | is a function of aerosol particle size, 
concentration and wavelength. The scattering efficiency 


factor weights the particle number differently for a given 


wavelength, so the radiance received at a given wavelength 


has’ information *eneethe mumbernsedistwaiaution ofmeepareiclesm 
This is the basis for the phase function parameterization 
used in the optical depth retrieval described in Chapter IV. 
Optical depth is a function of@@he extmmerionmeoe fates ene 
so it follows that changes in aerosol particle size or 
distribution will be reflected in optical depth and 


satellite detected radiance. 


B. RADIATIVE TRANSFER SOLUTION 


The general form of the radiative transfer equation for 
a given solar wavelength in a plane parallel atmosphere 


(arter Liou (P9ee))e as given by: 








dL, (0; wu, 9) O 
he Oil ee Lon Or ——_ || L, 0, Opie ae 
do AN ax 
oa) 
O, -¥, 
es ToD Ge ale. |e eee 
where uw is the cosine of satellite zenith angle (9), wo is 
the cosine of Sun zenith angle (90), » is the relative 


azimuth between the satellite and Sun, Lo is “ditrruce 


to 


radiance (W/m* um sr), Q is the solid angle defined by 9 and 
©, @ is the single scattering albedo, p is the scattering 
phase function, and F, is the incoming solar radiance at the 
top of the atmosphere. The total radiative energy in a beam, 
the energy scattered into the beam due to multiple 
scattering, and the energy scattered into the beam due to 
single scattering are the three terms on the right hand side 
of equation 2.4, respectively. 

Multiple scattering is considered to be small in 
atmospheres with small optical depths such as the cloudless, 


marine atmosphere. The single scattering approximation 


(after Liou (1980)), given by equation 2.5, results when the 


multiple scattering term is removed. 


eos oe ) ev.) = oH) ea 


an + p, 


tie aLgumentestot Ehe phase function (palo mine. scattering 
angle (w,). 
Durkee et al. (1991) show that in the cloud free, 


marine environment 6, is small enough to reduce Equation 2.5 


LO: 





ie = p(y, )S, (2.6) 


where the subscript “a” denotes aerosol related quantities. 
This linear relationship is obtained after removing 
contributions due to ozone, Rayleigh and ocean surface 
effects. This relationship also accounts for some multiple 
scattering especially where optical depth values are low 
(below 670.5). 


Aerosol optical depth (6,) can be obtained using 


a 
equation 2.6 if the remaining terms can be measured or 
calculated. Radiance at the satellite radiometer (L,) due 
to aerosol scattering at a given wavelength is measured. 
The incoming solar radiance (F,) 1s estimated given the sun- 
Farth geometry. The satellite zenith angle (represented by 
u) is determined by the satellite-Earth geometry. The 
scattering phase function (P) is parameterized (see Chapter 


iy) The single scatter albedo (@,) 1s the ratio of the 


oO 
scattering coefficient to the extinction coefficient. As 


Gieenesea earlier, the extinction CoerfIG¢Ient is the sum 


effect of absorption and scattering. Marine aerosols are 
weak absorbers therefore the scattering and extinction 
coefficients are approximately equal for this environment. 
This results in a single scatter albedo value of 
approximately one. However, there are some instances, such 
as when encountering combustion (carbon) and dust aerosols, 
when the single scatter albedo may be Significantly less 
than one. Aerosol radiance (L,) and scattering phase 


function (P) are described in detail below. 


adie Aerosol Radiance (L,) 


The radiance measured by a satellite radiometer results 
from a combination of sources including Rayleigh (molecular) 
scattering, surface and sub-surface reflectance, sun glint 
and aerosol scattering. To isolate the measured radiance 
due only to aerosol scattering the additional sources must 
be removed. The radiance sources can be accounted for in a 


linear expression (after Gordon and Clark, 1980): 
L=L- tL, - (1, + L,\t (2.7) 


where L, 1S radiance due to Rayleigh Scatter, L, is surface 


s 


metleectCed “radiance, Lr is (sun e@eint Sracdtane ey irc ees 


transmittance of the atmosphere. For satellite radiometer 
window channels in solar wavelengths, tT can be approximated 
as one with less than 10% error. 

A correction to L, is needed to account for ozone 
absorption in the upper atmosphere. L, 1s divided by the 
slant-path transmittance due to atmospheric ozone as in 


Equation 2.8: 


10 


i 
L = ———__ (2.8) 


1 a 
ses 
ai 
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where 6,, 1S ozone optical depth. Radiance contribution 
associated with Rayleigh scatter (L,) 1s calculated using a 
two-stream model by Turner (1973). Surface reflected 
radiance comes from several sources. Surface reflectivity 
over the ocean, however, is minimal since the ocean is 
assumed to behave like a blackbody. This assumption, 
peesentaod in Ramsey (4963), is based on ocean curiace 
reflectivity less than 5% for red-visible wavelengths. The 
MagenecOntributaen to Leis surface reflectance from foam 
(white caps) and sub-surface reflectance from suspended 
particles/plant life. Koepke (1984) describes foam 
reflectance as a function of wind speed, foam age and foam 


coverage. Ignatov et al. (1995), based on Koepke (1984), 


calculate foam reflectance (at 0.63 um) increasing from zero 


to 0.2% as the wind approaches 10 m/s. Additionally, 


Ignatov et al. (1995) calculate sub-surface reflectance (at 


0.63 um) in the open ocean, or coastal areas off the 
continental shelf in the absence of terrigenious influence, 
on the order of 0.14. This study approximates the combined 


effect of foam and subsurface reflectance as 0.5% at 0.63 um 
mage zerOolat 0.86 fm. Directional vet lectanee uemro hc 


variation in refractive index between the atmosphere and 
ocean is treated in the phase function term described below. 
€CentAtiw@ationgby spceular ug) 1 ni La omerhe scum iaicemen 
the ocean for low Sun angles can causes problems by 
artificially raising measured radiances. Measured radiances 
are not used if they are suspect based on the satellite- 
Barth-Sun geometry and wind speed. i cee Pree elec ae ee 


is assumed to be zero. 


Ia 


2. Scattering Phase Function (P) 


The scattering phase function determines in which 
direction the radiation is scattered. The scattering phase 
function is dependent on the radiation wavelength, but more 
importantly it is dependent on the size, composition, and 
distribution of atmospheric constituents. For Rayleigh 
scattering, the phase function is well understood. FOr 
particles whose circumference is near that of the incident 
radiation's wavelength, aerosols in this case, Mie theory is 
used to quantify the scattering phase function. Knowledge of 
the aerosol size distribution and composition is necessary 
to accurately quantify the phase function, but it is 


precisely the aerosol characteristics desired as a result of 


this retrieval method. To parameterize the phase function, 
several methods may be used. A common approach, described 
by Shettle and Fenn (1979), incorporates a phase function 


empirically fit using measured averaged size distributions 
and characteristics for a region. Another method, used by 


Durkee et al. (1991), takes advantage of the sensitivity of 


the phase function to radiative differences between two (or 
more) wavelengths to parameterize the phase function. This 
1s the method used in this study and is described in more 
detail in Chapter IV. 

There are three scattering paths that we are concerned 
with when using a single scatter model and they must be 
accounted for with the phase function. The pLrobaoll Tey er 
direct scatter off an aerosol particle back to the satellite 
sensor is determined by using the back-scatter angle (90- 
180°) as the argument for the phase function. The other two 
paths are considered diffuse scatter reflecting off the 
océan surface (bi-directional reflectance) . These include 


forward scatter (0-90°) off an aerosol particle toward the 


ocean and reflection back to the satellite sensor, and 


eZ 


reflection off the ocean surface and forward scatter off an 
aerosol particle toward the satellite sensor. The diffuse 
scatter can be estimated using Fresnal reflection 
coefficients. This leads to an equation for the effective 


phase function: 


Pert a a pale, a r, | (2 39) 


where P.,, 1s the effective phase function, P_ is the value 
of the phase function at the back scattering angle, P,is 
the value of the phase function at the forward scattering 
angle, and r is the Fresnel reflection coefficients at both 
Uu, and uw. Fresnel reflection coefficients are given in 


Bouatiem (2.10) : 


where 9; is the angle of incidence and 9, is the angle of 
transmission. From Snell's law, 9, = sin™’(sin(0,)/m)) where 


mis the index of refraction of seawater (1.33). 


dks 


ent 


Surface 


ome : . : oY 
SORT Cin en Oe EO ORE Rinne inh ee 





Figure 2.1. Radiative Transfer in the cloudless, 
marine environment. On the left, a typical 
atmospheric temperature profile depicting the marine 
atmospheric boundary layer (MABL) is included. Water 
vapor and aerosols are assumed to be confined to the 
MABL. Various scattering paths of incoming solar 
radiation (F,) from Sun to satellite are shown in the 
right panel. Radiance variations along these paths 
measured by satellite radiometers (L,) are primarily 
caused by upper atmosphere ozone absorption, Rayleigh 
(L,) and aerosol (L,) scattering, and ocean surface 
reflectance (L,) (Brown 1997). 
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IIrt. DATA 


In order to validate the optical depth retrieval method 
described in this study, reference data from ACE-2 were 
chosen based on the availability of surface and airborne 
sunphotometer measurements of optical depth. PO -aaiehie 
greatest extent possible, sunphotometer measurements are 
matched to satellite observations both spatially and 
temporally. The aerosol optical depth composites were 
produced using the validated retrieval method with AVHRR 
data collected during each of the three experiments. This 
chapter will briefly describe the data sets and 


instrumentation used to collect the data used in this study. 


A. OVERVIEW OF EXPERIMENTS 

The main goal of ACE-1 was to determine and understand 
the properties and controlling factors of aerosols in the 
remote marine atmosphere that are relevant to radiative 
honcingsandesclimate:. This experiment took place in the 
minimally polluted marine atmosphere of the southern ocean 
south of Australia from November 15 to December 15, 1995. 
This area was selected due to the relatively simple marine 
aerosol background and its distance from the Northern 
Hemisphere sulfate aerosols. This aerosol background can 
provide a baseline to compare with anthropogenically 
perturbed aerosols transported from the continents. Data 
was collected at two land based sites (Cape Grim, Tasmania 
and Macquarie Island), two ships (Australian R/V Southern 
Surveyor and NOAA R/V Discoverer), and one airplane (NCAR C- 
ino 0) - Full resolution NOAA-14 AVHRR data was downloaded 


iS 


real-time in Hobart. Bates et al. (1998) provides an 
overview of ACE-1. 

TARFOX was designed as a closure study to better 
understand the radiative forcing effects of aerosols. 
TARFOX was conducted in the continentally influenced 
environment off the eastern coast of the United States near 
Waltlops, Island, Virginia from 10-31, Juig@ 1996. During 
TARFOX, a variety of aerosol conditions ranging from 
relatively clean to moderately polluted were observed and 
measured. In situ measurements were conducted by airborne 
platforms including the University of Washington's (UW) 
C-131A and the United Kingdom's (UK) Meteorological Research 
Flight C-130 (see Russell et al., 1996). Satellite imagery 
collected included full resolution NOAA 14 AVHRR data. The 
TARFOX Operations Summary (Whiting et al., 1996) contains 
details of the field collectiony c#rorc. Brown (1997) 
provides a validation of the methods used in this study from 
TARFOX field observations. 

ACE-2Z “Carriedijonm the goals of @AeGE-1, but Bteeused on 
anthropogenic aerosols from the European continent and 
desert dust from the African continent as they move over the 
North Atlanetemeecean- This experiment was carried out 
between 16 June and 25 July, 1997 and involved coordinated 
data CoOlLe@CTION “Diese Sik QbuGrart, one, Ship, “and cqreund 
stations on Tenerife, Portugal and Madeira. An overview on 
the ACE-2 objectives can be found in the experiment’s 
Science and Implementation Plan (IGAC, 1995). NOAA-14 AVHRR 
data was collected in real-time during the exercise. 
Aerosol optical depth measurements from four different 
sunphotometers were used to validate aerosol optical depth 
values from the satellite retrieval. The NASA Ames Airborne 


Tracking Sunphotometer (AATS-14) was flown aboard the Center 


nC 


for Interdisciplinary Remotely Piloted Aircraft Studies 
(CIRPAS) Pelican aircraft. The NASA Ames Airborne Tracking 
Sunphotometer (an TS— 6) was operated aboard the R/V 
Vodyanitskiy. The ground station at Tenerife used a CIMEL 
Electronique 318A sun/sky scanning spectral radiometer and 
meine -GrOUNd Staten at Sagres, Portmgal usedmboerh the UVISIR- 
1 and IR-RAD sunphotometers. 


B. INSTRUMENTS 


ska NOAA Advanced Very High Resolution Radiometer 
(AVHRR) 


The AVHRR instrument iS a component of the NOAA Polar 
Orbiting Operational Environmental Satellite (POES) series 
satellites. These satellites are in Sun synchronous orbit 
(883 km) and provide two passes per day in the morning and 
evening, respectively. The NOAA-14 data were used 
exclusively for this study due to sun-earth-satellite 
geometry necessary for the aerosol optical depth retrieval. 
The AVHRR instrument measures radiant and solar-reflected 
energy from sampled areas of the Earth in five spectral 
bands with a sub-satellite resolution of 1.1 km. Table 3.1 
lists the characteristics of the individual radiometer 
channels for the AVHHR (Kidwell, 1995). Channels 1 and 2 
are used in the optical depth retrieval. Channels 1 through 
5 are used in the cloud screening analysis. 

All AVHRR channels are calibrated prior to launch. 
Channels 1 and 2 of the AVHRR have no onboard calibration 
systems. Post calibration methods for these channels have 
been developed by the NOAA/NESDIS Office of Research 
Applications based on the results of Rao and Chen (1995). 
Calibration formulae are incorporated into the satellite 


image processing discussed in Chapter IV. 


iy 
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Table 3.1. NOAA AVHRR Radiometric 
Channels 


2s NASA Ames Airborne Tracking Sunphotometer (AATS) 


By tracking the Sun, the NASA AMES Airborne Tracking 
Sunphotometer measures the relative intensity of the direct 
incoming solar radiation in multiple spectral channelise 
This information is converted to optical depth (Livingston 
and Russell, 1997; Schmid et al., 1999). The AATS-14 is the 
14-channel variant and.was flown aboard the CIRPAS pelican 
a iggC Gat tC Channel filters are at wavelengths from 380 to 
1558 "nme "By flying amwalrcraft at “Yow alritudes in ”cloud 
free regions within the atmospheric boundary layer, total 
column optical depth can be well approximated. FOr: Chis 
data set the Pelican flew mostly off the Island of Tenerife 
in cloud free areas. The AATS-6 is the 6-channel variant 
and it was operated aboard the research vessel R/V 
Vodyanitskiy. Channel filters are at wavelengths from 380 
op 20) alae The R/V Vodyanitskiy sailed mostly south of 
Portugalyemear the Str. of SGibralitar veelicceting data ane 
area influenced by anthopogenic aerosols. NASA Ames 
performed instrument calibration prior to operation. These 


instruments are designed to maintain calibration within 1% 


ine 


dum@ing Cperation® Resolution of optical depth” is on the 


order of 0.01. (Matsumoto et al., 1987). 


3. CIMEL Electonique 318A Spectral Radiometer 


The ABRONET program (see Holben, 1998) uses the CIMEL 
Electronique 318A sun/sky scanning spectral radiometer. 
This solar powered all weather instrument is robotically 
pointed and has approximately a 1.2 degree full angle field 
of view. The radiometer makes only two basic measurements, 
either direct sun or sky. The direct sun measurements are 
made in eight spectral bands between 340 and 1020 nm. The 
instrument was located at 28.0° N 16.6° W and an altitude of 
coe mon Tenerife is lane The network of instruments is 
Galibrated twice per year and obtains aerosol optical 
Moemness accuracy to less than 0.01 for waveMengtns Greate= 


jepemeeaa( nm. 


4. UVISIR-1 Sunphotometer and IR-RAD Sun-radiometer 

The UVISIR-1 and IR-RAD instruments were located.in 
Sagres, Portugal at 8.9° N 36.9° W and an altitude of 50 m. 
These instruments were operated by the Institute of Physics 
and Chemistry of the Upper and Lower Atmosphere (FISBAT), 
Bologna, Italy. The UVISIR-1 instrument scans at 12 
wavelengths between 360 and 1050 nm with a scanning interval 
Cu2 min. The IR-RAD instrument scans at 13 wavelengths 
between 400 and 3700 nm with a scanning interval of 3 min. 
peeneck or. the calibration constame reliability i1seean 
‘progress, therefore the aerosol optical thickness accuracy 
is estimated at 0.01 in line with the other sunphotometers 


iseceumn Chis study. 


ro 





IV. OPTICAL DEPTH RETREIVAL AND COMPOSITING PROCEDURES 


This chapter describes the procedures used to calculate 
aerosol optical depths from the satellite imagery and the 
procedures used to create the composite images for the three 


exercise areas. 


A. SATELLITE IMAGE RETRIEVAL/DISPLAY 


The NOAA POES High Resolution Picture Transmission 
(HRPT) data collected during the exercise periods were 
archived to 4mm data cartridges. The Terascan Earth Remote 
Sensing System (by SeaSpace Corporation) was used to 
retrieve and display this data at full resolution. Terascan 
system includes software that converts sensor radiance 
eounts Into albedo/brightness temperature and calculates 
required angles from Sun-Earth-satellite geometry A full 
suite of post processing enhancements is also included. 


This greatly aided in the optical depth analysis. 


B. OPTICAL DEPTH RETRIEVAL 


The optical depth retrieval technique used in this 
study 1S an automated process based on work done by Brown 
(1997). Terascan and FORTRAN 77 code are used to manipulate 
the satellite data and calculate the required information. 
Terascan commands are used to retrieve raw satellite images 
from the pass disk of the Terascan receiver system. 

The FORTRAN code included a check for sun glint, a 
cloud screening algorithm, and all the radiative transfer 
calculations including optical depth. The sun glint check 
determined a probability of sun glint in a given pixel based 
on the Cox and Munk (1954) model. This process takes into 
aceCeuns sun-satellite geometry and windé= speed. A 
conservative 14 m/s was used for the wind speed over the 


entire image and the pixel was removed if the probability of 
Zak 


sun glint was greater than 35%. mnie ™"c lemd "Screening 
algorithm, based on Saunders and Kriebel (1988), applies the 
following tests to determine cloud contamination in a given 
Diels: 

-Gross cloud check: The pixel was removed if the 
Channel 4 temperature was less than the sea _ surface 
temperature (additionally, if the channel 4 temperature was 
greater than 303 K the pixel was considered to be land). 

-Spatial coherence: The pixel was removed if the 
Standard deviation of the channel 4 temperature was greater 
than 0.1 K or the standard deviation of the channel 3 
temperature was greater than 0.45. 

-Dynamic reflectance threshold test: The pixel was 
removed if the channel 2 reflectance was greater than 15% 
and if the channel 4 minus channel 5 temperature difference 
is greater than zero. The temperature difference test allows 
for a high reflectance due to dust aerosols (used only for 
ACE-2). 

-VIS/NIR ratio test: The pixel was removed if the ratio 
6f Channel 1 tem channel! 2 bi-directional reflectance: was 
fess than 1,33) 

-Thin cirrus test: The pixel was removed if the 
difference between the channel 4 and channel 5 temperatures 
was greater than pre-computed clear sky values. 

The radiative transfer calculations are discussed 
below. After these calculations, the radiative transfer 
output from the FORTRAN code was imported back into Terascan 


format for post-processing and enhancement. 


Ge RADIATIVE TRANSFER CODE 


The radiative transfer solution (Equation 2.6) is 
calculated using FORTRAN code and follows the theory in 
Chapter II. Optical depth is the primary output along with 
the solution to Equation 4.2. The input wavelength used is 


the central value of the satellite channel's spectral band. 
Ze 


ie Calculation of specifiegwinput variables “for "quataen 2.6 


is described below. 


Le Solar Radiance, Ozone and Rayleigh Optical Depths 


Input values for solar irradiance are taken from the 


NOAA Polar Orbiter Data Users Guide (Kidwell, 1995). The 


values are determined by calculating a weighted average of 
solar irradiance over the spectral band of satellite 
radiometer. Solar radiance is calculated by dividing solar 
to@agiance, by piwis./1) . E, is corrected for variation in 
the Earth-Sun distance prior to input into the model. Input 
values for ozone and Rayleigh optical depths are based on 
those reported Elterman (1970). Table 4.1 lists the values 


of solar irradiance, and both ozone and Rayleigh optical 


depths used in these calculations. 


ee 


NOAA 14 AVHRR Ch 1 C027 O87 LeZ2e 


NOAA 14 AVHRR Ch 2 OG O oi 1030 


Table 4.1. Values of Ozone and Rayleigh Optical Depths, and 






Solar Irradiance. 


Ze. Scattering Phase Function 


Mie theory can be used to calculate phase function if 
the aerosol size distribution is known. The aerosol 
distribution is not normally available unless in situ 
aerosol measurements are taken. Parameterization of the 
scattering phase function is necessary to solve _ the 
radiative transfer problem. The method used in this study 
takes advantage of the differences in the measured radiance 
at different wavelengths to parameterize the scattering 


phase function. 


OM 


Durkee et al. (1991) proposed using the measured 


aerosol radiance differences in the AVHRR channel 1 
(visible) and 2 (NIR) to parameterize the phase function. 
The scattering efficiency (Q,.,,.) of an aerosol distribution 
is wavelength dependent, therefor, when the radius of the 
aerosol is nearly equal to the radiation wavelength, Q,.,, is 


a maximum. The ratio of the channel radiances will be larger — 


for smaller size particle distributions and smaller for 


larger size particle distributions. Durkee et al. (1991) 
called this ratio the particle size parameter, S,.. Since 
S,, varies in each pixel of the satellite image, the 


scattering phase functions can be parameterized for each 
pixel, allowing variations in aerosol distributions in the 
optical depth retrieval. 

Seven model aerosol size distributions (M0O-M6) 
(based on Brown, 1997) were created. The scattering phase 
functions and extinctions for these models were calculated 
using Mie theory. The models include one single-mode log- 
normal omic! gull evens oval and S1x bi-modal log-normal 


distributions which follow the equation: 





dN (r) N; exp | 1_( r= in rs 
woe = me — SS Ao 
ee SiC 207\  in(10) a 
where j} = 1,2, N is the number density, r, is the mean 


radius and o is the standard deviation of log(r). 


The mode radii and standard deviations used are 
designed to model the typical variations of aerosol 
dTSEmiouelens  lasric Nabine aanespoi- Goi First mode im 
the distribution models of the background aerosol while the 
second mode in the distribution models ocean-produced 
aerosol. Table 4.2 lists the mode radii and standard 


deviations. Figures 4.1, 4.2, and 4.3 illustrate the 
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resulting aerosol size distributions and the corresponding 


scattering phase functions at 0.63 um and 0.86 umm. 


Model | Mode Radii (pm) Number Std Dev (o) 
Density (N) 
















Table 4.2. Mode Radii, Number Density and Standard 


Deviations for Model Aerosol Size Distributions 
Si. values are calculated using the model phase 


functions and extinctions in the following equation (based 


on linear single scattering theory) : 


Hi fe eee eeu 
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(4.2) 
Lone Reno oa © extenz 


S,, 1s corrected for water vapor absorption. in AVHRR channel 
2 after Mahony (1991) using the split-channel (channels 4 
and 5) water vapor retrieval proposed by Dalu (1986). # The 
resulting S,, values for each aerosol model are shown in 
Figure 4.4. 

Si. values are calculated for each pixel from AVHRR 
channel 1 and channel 2 data. Computed S,, values along 
with scattering angle are entered into a look up table (LUT) 
that represents Figure 4.4. The model aerosol distribution 
(Figure 4.1) best represented by the observed radiance 


measurements is selected. By using the scattering angle and 


a 


selected model size distribution, phase function values are 
then selected from LUTs based on Figures 4.2 and 4.3. 

An additional check was added for extremely clean 
environments. In some cases during ACE-1 and ACE-2 the 
measured aerosol radiance was so low that it was almost 
indistinguishable above the background (or bias) level. 
This resulted in an artificially high S,, and an incorrect 
aerosol model assumption. In these cases, defined by 
brightness counts less than 2.5 above bias, the aerosol 
model was set to M6. The M6 model is more typical of this 
clean environment. Optical depths were than recalculated 


based on the new phase function selection. 


D. COMPOSITE PROCESS 


The composite process consisted of grouping pixel values 
iiniie. 51 Nnseandwpe sf omnes bataciaca 1. calieudlatoms sees Ooi 
km bins were defined across a predetermined area of 
coverage. These areas were slightly larger than the 
exercise areas themselves. The values calculated by the 
radiative transfer code for each NOAA 14 pass during the 
exercise were summed and averaged over a bin. In _ some 
instances a Standard deviation was calculated. Since the 
AVHRR resolution iS approximately 1 km, there could be as 
many as 100 pixels per bin. Some manual and automated 
quality control was used. Fach pass was manually screened 
for obvious errors such as optical depth values within the 
boundaries of a cloud field. One hundred pixels from each 
side of the image are removed to avoid sampling data from 
the edges of the pass where distortion exists. Summing all 


of the passes then created the area composite. 


Zo 


dN/dr 


Model Aerosol Size Distributions 
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Radius (um) 
Model Aerosol Size Distributions. 
lognormal aerosol 


distributions with MO representing the background 


Figure 4.1. MO-M6 
correspond to model two-mode, 


aerosol mode and M6 representing the largest 


(continental) 
(From Brown 1997). 


oceanic aerosol mode 
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V. RESULTS 


A. VALIDATION 


The satellite retrieval validation was accomplished by 
comparing the retrieval aerosol optical depth (AOD) output 
to sunphotometer AOD data collected during ACE-2. Five 
different instruments collected sunphotometer data at 
various times and places throughout the exercise period. 
The sunphotometer data used was correlated both spatially 


and temporally with the retrieved data to the greatest 


extent possible. Comparisons were made at 630 nm and 860 
nm. In all cases the sunphotometer data was linearly 
interpolated to these wavelenghts. Only the NOAA-14 


afternoon passes provided the proper sun-satellite geometry 
necessary to calculate AOD using this retrieval method. The 
uncertainty associated with the resolution of the satellite 
radiance measurements results in an average AOD uncertainty 
of + 0.02 based on Brown (1997). ihe Yesul eanibeceita yee 
can be seen in Table 5.1. 

A CMEL Electronique 318A sunphotometer, operated by 
AERONET, was located at the Tenerife ground station. This 
station was located near Punta del Hidalgo at 28° N 16.6° W 
Boge webevatileon Of ~L0  m- Retrieved AOD data from the 
nearest pixel just off the coast from this ground station 
was used in the comparison. The sunphotometer data set was 
prescreened for clouds by AERONET. Bight matches to 
retrieved AOD data were found between 30 June and 18 July. 
The time difference between retrieved values and 
sunphotometer values ranged between zero and 39 min. An 


error of + 0.01 is associated with the AOD measured by this 


instrument (Holben 1998). 


on 
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a0-sanl 150] 0-08] 0-08 
a=suil 2509]_0-0e| 0.05] 180[0,os6[0-037[FENERTFE [0000 
[25-Sun| 1434] 0.1 0.08] 1434/0. 1270. O78|AATSE 
0.125 
[22-Fun| 1507] 0.08] 0.06| 1508|0.078|0.053|IRRAD- SAGRES [0001 
[23-Sun| 1456|0.076|0.066| 1454] 0.066 0.06|IRRAD SAGRES [0002 
[24-Tun| 1445 0.06 0.05] 1447[0.055[0.037/TRRAD SAGRES [0002 


aS 5 = Valndaeion ae Bee sr the a 
retrieval data and sp is the soso eae ei 


The ground station at Sagres, Portugal was located 


at the coastal military base Radio Naval de Sagres, 37° N 


8.9° W and elevation 50 m. Two sunphotometers, UVISIR-1 and 
IR-RAD, were operated at this station by FISBAT. Retrieved 
AOD data from the nearest pixel just off the coast from this 
ground station was used in the comparison. Two matches to 
retrieved AOD were found for the 11 and 12 July data from 
the UVISIR-1 instrument. Seven matches to the retrieved AOD 
data were found between 22 June and 12 July from the IR-RAD 
instrument. Errors for the data collected by the two 
sunphotometers at this station were not available, but were 


assumed to be similar to the error reported for the other 
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instruments used during the exercise and assigned a value of 
aoe OI 
NASA AmeS operated the  AATS-6 aboard the R/V 


Vodyanitsky. The ship spent most of its underway time south 
of Sagres outside the Strait of Gibraltar, but in one 
coordinated experiment the ship traveled farther south off 
tiemcOast Of Morocco: AOD values from this sunphotometer 
were screened for clouds and ship mast interference. Two 
matches were found, one on 24 June and one on 25 June, 
during the time when the ship was outside the Strait of 


Gibraltar (35.6° N 9.2° W and 36.1° N 9.0° W respectively). 
The temporal agreement in these cases is good. One match 
was found for 10 July when the ship was off the coast of 
Morocco (29.3° N 11.9° W). There is a 21 min. difference in 
the data times for this match, but this is considered 
acceptable since the ship moves a very short distance and 
the aerosol loading can be assumed to be consistent over 
this time. The error associated with the AATS-6 AOD values 


are between + 0.006 and + 0.008. 


The AATS-14 was flown aboard the CIRPAS Pelican 
aircraft. The majority of flights were in clear air masses 
nearby the Island of Tenerife, but one flight was near the 
coast of Morocco in order to coordinate with ship 
measurements. The aircraft flew at many altitudes due to the 
requirement of other instrumentation onboard. The AOD 
values used for comparison were only those calculated when 
the plane was at an altitude of 100 m or less and near the 
time of the satellite retrieval. For this short period of 
time (6 min or less) the AOD values were averaged and the 
average was used as a Single value for comparison. Three 


Mamehnes were tfoune@®s for this data set:21 Jmne meare2zs.9° N 
17.2° W had a temporal difference of less than 11 min, 10 


July near 29° N 11.13° W had a temporal difference of less 
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than 13 min, and 17 July near 27.8° N 16.6° W had a temporal 
difference of less than 34 min. The error associated with 
the AATS-6 AOD values are between + 0.004 and + 0.008. 


Some retrieved AOD data was not used even though there 
was matches to sunphotometer data. In areas contaminated by 
sun glint or clouds, retrieved AOD values were automatically 
throw away by the retrieval process. AOD values that were 
near the edge of the satellite pass (defined by 100 pixels 
in the composite process) were not considered valid. 
Finally, retrieved values that had valid sunphotometer 
matches but were in a strong gradient of AOD were not used. 
Strong gradients of retrieved AOD, doubling or tripling of 
the AOD value over a few pixels, are suspect for matches 
because they may be inconsistent both spatially and 
temporally with the sunphotometer AOD and it is difficult to 
place the match within the gradient. 

The scattering angle used in the calculation of 
HebwaleVvedmnOl was Checkecdmro ensue |Eheue was nomoras stoma 
single angle when choosing a value for the phase function. 


As seen in figures 4.2 and 4.3 the value of the phase 


function converges near a scattering angle of 140°. Data 
with scattering angles near 140° are not affected by aerosol 


model type and would therefore not provide a robust 
validation. The range of scattering angles for this data 


Seneis betweem e152 “and 9276 Swhich 1S Gn anweacoaaweon 


Significant spreading of modeled phase function values. 

The resultant data set is shown in a scatter plot in 
Figure 5.1 and indicates a positive validation of the 
retrieval method. Sunphotometer source and wavelength data 
points are grouped together in the plot. There is very good 
agreement in the lower optical depth range for all sources 
and both wavelengths. In the higher optical depth range 


the correlation is slightly off, more so for the values at 
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860 nm. Since these high AOD values are primarily due to 
dust aerosol, the bias may be due to the non-absorbing 
aerosol assumption made in the satellite retrieval method. 
If aerosol absorption was included the retrieved AOD values 
would be higher and correlate more closely with the 
sunphotometer AOD. Dust aerosols typically have single 


scatter albedo values of 0.9 or less (Ignatov et al., 1995). 


Even with this problem the correlation coefficient has a 
value of 0.88. Additionally, the standard error for the 630 
nm and 860 nm data set is 0.023 and 0.022 respectively. 
Brown (1997) using data from TARFOX, also achieved a 
positive validation of this AOD retrieval method. ‘The 
validawion of the AOD retrieval method»1is assumed, for athe 
ACE-1 composites discussed below. 

Another validation method used by Nakajima and 
Higurashi (1998) compares the Angstrém exponent calculated 
from satellite retrieved AOD with those calculated from 
sunphotometer AOD. A scatter plot of the Angstrém exponents 
calculated from the ACE-2 data set 1S given in Figure 5.2. 
There 1s good agreement between the calculated values with 
the exception of a few outliers. The outliers are all 
points from the Sagres data set and all had extremely low 
values of AOD at both 630 nm and 860 nm. Optical depth 
values this low are calculated from radiances near the 
detection limit for AVHRR where the magnitude of AOD is 
about the same as the uncertainty due to digitization of the 
radiance. This can lead to anomalously high or low ratios 
of AOD and therefore anomalous values of the Angstrém 
exponent. The two outliers in the upper left of figure 5.2 
are both from 11 July and both have high satellite retrieved 
AOD ratios compared to the sunphotometer ratios (see Table 
5.1 for data values). The third outlier, on the right side 
of Figure 5.2, is from 30 June. The satellite retrieved AOD 
are slightly higher than the sunphotometer AOD, but their 
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relative difference is similar. Since the AOD values are low 
the ratio of the sunphotometer AOD values produces a higher 
Angstrém exponent. The standard error of this data set with 
the outliers is 0.73 and without the outliers is 0.41. This 


adds to the confidence of our validation. 


B. REGIONAL ANALYSIS 


The regional composite images can been seen in figures 
5.3-5.7. The composites cover larger areas than the actual 
exercise areas in order to identify possible sources and 
transport of aerosols. Figures 5.3 and 5.4 (a)-(c) are 
images of AOD, at 630 nm and 860 nm respectively, for the 
three exercise areas. Figures 5.5(a)-(c) are composite 
images of AOD ratio. The ratio of AOD at two different 
wavelengths is related to the Angstrém exponent and is an 
indicator of aerosol size distribution. Higher Angstrom 
exponent values indicate areas of pollution or biomass 
burning and lower values indicate marine and dust aerosols 
(Nakajima and Higurashi, 1998). Figures 5.6(a)-(c) are 
images of the standard deviation for AOD at 630 nm. These 
images are used to indicate the consistency of the aerosol 
properties over the exercise period. The number of pixels 
per bin that were used in the composites is shown in figures 
Sey (avec). Higher numbers correlate to areas which @wheme 
generally cloud free and not affected by sun glint. Numbers 
above 100 represent a contribution a more than one satellite 
pass used in the composite. The statistics cited below are 
for the exercise areas only. A white box on the images in 


figures 5.3(a)-(c) indicate the exercise areas. 


us ACE-1 


The ACE-1 study area, defined by research aircraft and 


ship operations, was between 40°-55°S and 135°-160°E. During 
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the field operations period there was an above average 
occurrence of cold frontal passage. A long wave trough 
passed over the exercise area during the course of the 
experiment. This led to strong westerly to northwesterly 
flow in November and southerly to southwesterly flow in: 
December. Sea salts dominate aerosol properties in this area 


(Hainsworth et al., 1998). As expected the AOD calculated 


(figures 5.3(a) and 5.4(a)) were low over most of the area 
with a mean of 0.13 at 630 nm and 0.11 at 860 nm. This is 
consistent with the assumption that continental aerosols are 
not the primary influence in this area. There are some 
areas of higher AOD over the southern ocean near the sub- 
tropical convergence Zone. This may be due to increased 
wave action and therefore higher surface reflection 
contamination in the retrieval. Additionally, this area is 
also known to be a large producer of non-sea salt sulfates 
that contribute to aerosol concentrations and may lead to 
higher AOD. Individual images of retrieved AOD from early 
in the exercise period revealed higher AOD off the East 


Coast of Australia consistent with continental aerosols 


being advected by the observed westerly winds. Some high 
AOD ratio values (figure 5.5(a)) are scattered over the 
eastern portion of the composite. This may indicate an 


influx of anthropogenic aerosols into this area or an 
abundance of biologically produced sulfate aerosols. The 
Standard deviation of retrieved AOD (figure 5.6(a)) is low 
with a mean of 0.028 at 630 nm and 0.064 at 860 nm. This is 
indicative of the relatively homogenous aerosol properties 


in this area. 


Ze TARFOX 


TARFOX was conducted in the continentally influenced 


environment off the East Coast of the United States near 


a 


Wallops Island, Virginia. A variety of aerosol conditions, 
from relatively clean to moderately polluted, were observed. 
A persistent upper-level trough created extensive cloudiness 
on many days and highly variable haze conditions due to 
frequent shortwave weather systems passing over the area. 
Additionally, Hurricane Bertha threatened the exercise area 
for three days. The TARFOX composites extend south and west 
in order to identify aerosol sources from sub-tropical 


latitudes and aerosol transport across the North Atlantic 


Ocean. AOD composites are shown in figures 5.3(b) and 
5.4(b). The mean values from these composites are 0.353 at 
630 nm ANG. 2 7deeet eee nm. These higher values are 


consistent with continental soil and pollution aerosols. 
Au@ener se courcemlcmeminerd lyedust = that Jorigqinate cme) nesltosee 
Africa. This dust makes it way across the sub-tropical 


North Atlantic Ocean to the Caribbean where it gets 


entrained in the mid-latitude westerly flow. This process 
1s evident in the composite images. Extensive cloudiness 
may bias the retrieved AOD. The cloud screen algorithm has 


difficulty near cloud boundaries and it is in these areas 
that aerosol humidification may cause higher calculated. 
valueseet AGD. gbhe AGD ratiao..(figure 5.5(b).as highs near 
areas where pollution aerosols are expected, specifically 
off the coast of major cities in the North East U. S. This 
feature, due mainly to the increased amount of small aerosol 
particles from combustion processes, persists well into the 
North Atlantic Ocean. The mean values of AOD standard 
deviation for this region were 0.148 at 630 nm and 0.112 at 
860 nm. The relatively high standard deviation of AOD 
(figure 5.6(b)) 1s indicative of the variability of the 


aerosol conditions associated with this region. 
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a ACE-2 


The ACE-2 area was selected based on its exposure to 
various aerosol types, including clean marine, anthropogenic 


from Europe and desert dust. The main study area was 
between 23°-44°N and 8°-25°W. The composite area extends 


further west into the North Atlantic Ocean and east into the 
Western Mediterranean Sea. The mean AOD values are 0.142 at 
630 nm and 0.126 at 860 nm. This is consistent with most of 
the area being dominated by clean marine aerosol. The most 
Significant feature is the long plume of high AOD that 
extend off the African continent into the North Atlantic 
Ocean (See figures 5.3(c) and 5.4(c)). AOD values 
approaching 1.0 are common in this plume. An anomalous area 
of high AOD is present in the northwest corner of the 
composite. This may be an extension of the dust plume that 
has been caught in the mid-latitude westerlies or aerosols 
that made there way from the Western Atlantic around the 
sub-tropical high. AOD Plumes extending off the coast of 
Algeria and Tunisia into the Mediterranean Sea _ are 
associated with dust carried by Sirocco winds. There is 
some evidence of pollution out breaks, based on higher AOD, 
south and west of Portugal. These events are seen 
periodically in individual cases, but are averaged out over 
time in the composite. The mean standard deviation of AOD 
was 0.064 at 630 nm and 0.061 at 860 nm. The AOD ratio 
(figure 5.5 (c)) is low and fairly consistent over the 
entire area. North of the dust plume the ratio values 
increase which may be a result of pollution aerosols from 
the European continent. The mean standard deviations are 
biased by large standard deviation values in the area of the 
African dust plume (see figure 5.6(c)). 

A scatter plot of Angstrém exponent versus AOD is 


useful for an optical characterization of aerosols in terms 


oy 


of the typical particle size and distribution. @@#he typical 
observed value of the Angstrém exponent is approximately 1 
for small urban-type aerosols and approximately 0O for large 
soil-derived particles (Nakajima and Higurashi, 1998). 
Figmiees Sacmes Sa scatter plot of Angstrém exponent vs. AOD -: 
for the ACE-2 data set. The Sagres and AATS-6 data points 
tend to have higher Angstrém exponents, although lower AOD, 
as is expected due to urban-aerosol effects from the 
European continent. The AERONET and AATS-14 data points 
tend to have lower Angstrém exponents, with varying AOD, as 
is expected in the open ocean near the African dust plume. 
The satellite data points correlate well with the respective 


sunphotometer data points. 


4. Regional Comparison 


Regional comparison is accomplished through frequency 
plots of AOD at 630 nm and 860 nm wavelenghts (figure 5.9). 
The data used in the histogram is limited to only the 
operational areas of each exercise as described above. The 
histogram plots for ACE-1 and ACE-2 have a similar profile 
with modes of 0.115 and 0.095 (at 630 nm) respectively. 
This is indicative of the relatively clean environment in 
these regions. The effect of the African dust plume in ACE- 
2 can be identified by the tail of the histogram plot toward 
higher values of AOD. It is interesting that the ACE-2 
frequency plot has a lower mode than the ACE-1 plot even 
though the mean AOD in the ACE-1 composite is larger. This 
1s a good example of regional differences in aerosol 
distributions. The ACE-1 area is influenced by essentially 
one type of aerosol constituent that leads to a similar mode 
and mean. The ACE-2 area has Several aerosol constituents. 
The main constituent at the dominant mode is similar to ACE- 
1, but its lower mode value implies smaller particles. The 


other constituents, dust and anthropogenic aerosols, 


40 


contaminate the area and lead to the higher average AOD. 
The TARFOX plot has a much higher mode (0.185 at 630 nm) and 
a Significantly larger spread. This is indicative of the 
larger aerosol loading and varied aerosol types in this 


region. 
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Satellite Optical Depth Retrieval 


Figure 5.1. 
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Surface Optical Depth Measurements 
Scatter diagram of surface measured aerosol 


optical depth versus aerosol optical depths retrieved from 


satellite. 


The correlation 1s very strong and the standard 


error is low. Absorption by aerosols is considered to be 
zero in the retrieval method. This may account for the 
weaker correlation at higher values of optical depth. 
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Qa satellite 
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0 1 2 3 4 


a sunphotometer 
Figure 5.2. Scatter diagram of Angstrém exponents (a) 
calculated from sunphotometer AOD vs satellite retrieved AOD. 
There is good agreement between the two data sets for low 
values of a. Small values of @ indicate an area of soil- 
derived particles such as the ocean area next to the Sahara 
desert (Nakajima and Higurashi, 1998). The few outliers that 
are present are due to low values of AOD at these points 
which leads to a higher error. The standard error for this 
data is 0.73 with the outliers and 0.41 without the outliers. 
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Figure 5.9. Histograms of aerosol optical depth at 630 nm and 
860 nm for ACE-1, TARFOX and ACE-2. 
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VI. CONCLUSIONS/RECOMMENDATIONS 


re CONCLUSIONS 


A radiative transfer algorithm in the solar wavelengths 
for the NOAA POES AVHRR is presented for the cloud-free, 
marine atmosphere. Priore COnteimneOMomdies linearized, 
Single-scattering theory with an estimate of bi-directional 
surface reflectance. Phase functions are parameterized 
using an aerosol distribution model and the ratio of 
radiance values measured in channels 1 and 2 of the AVHRR. 
The aerosol distribution model is based on a _ bimodal 
lognormal distribution with significant modes at 0.1 and 0.3 


lim. Ratios o£ channel 1 to channel 2 are crossed with the 


scattering angle in a look-up table (LUT) to determine the 
best representative aerosol distribution. Using this model 
and, again, the scattering angle the phase function is 
chosen from a LUT calculated from Mie theory. Automated 
cloud screening iS accomplished through preset thresholds 
and known comparisons using AVHRR channels 2-5 based’ on 
Saunders and Kriebel, 1988. Sun glint removal is based on 
Cox and Munk (1954) and takes into account sun-satellite 
geometry and wind speed. 

Retrieved satellite aerosol optical depth (AOD) was 
compared to surface measured sunphotometer AOD collected 
during the International Global Atmospheric Chemistry (IGAC) 
Project's Second Aerosol Characterization Experiment (ACE-2) 
from June 16 to July 25, 1997. Sunphotometer data from four 
independent sites were used in the comparison, ground 
stations at Tenerife Island and Sagres, Portugal, AATS-6 on 
board the R/V Vodyanitskiy, and AATS-14 on board the CIRPAS 
alrcraft Pelican. There was very good agreement between the 


AOD retrieved from satellite and those measured by 


ol 


sunphotometer, especially at low values of AOD. At higher 
values of AOD the retrieved values are low, but this is most 
probably due to the no-absorption assumption made in the 
retrieval. The comparison data set has a correlation 
coefficient of 0.88 with a standard error of 0.02 at both 
channel 1 and 2 wavelengths. This indicates a positive 
validation of the retrieval method. Another validation 
method, used by Nakajima and Higurashi, 1998, compares 
Angstrém exponents calculated from the data set. A similar 
comparison done here shows good agreement between the 
retrieval AOD and the surface measured AOD. This adds to 
the positive validation. Validation of this retrieval 
technique was also seen in Brown (1997) using data from 
TARFOX. The retrieval technique is assumed to validate for 
the ACE-1 area and is used for the regional analysis 
discussed below. 

Regional aerosol properties are examined with an 
emphasis on the differences between the ACE-1, TARFOX and 
ACE-2 regions. AOD determined from satellite retrieval were 
averaged for the duration of the exercise periods and 
combined into composite images at 10X10 km resolution. 
Statistics were also calculated for each exercise area. The 
images give a good indication of aerosol loading as well as 
source and transport regions in each area. A scatter plot 
of AOD versus Angstrom exponent for the ACE-2 data set is 
used to examine aerosol characteristics. Data from Sagres 
and R/V Vodyanitskiy have higher Angstrém exponents and 
lower AOD consistent with urban-type aerosols. Data from 
Tenerife and Pelican have lower Angstrém exponents and a 
varying AOD consistent with dust aerosol influence. AOD 
frequency plots are used to identify region differences of 
the three exercises. ACE-1 and ACE-2 regions have strong 
modes at AOD at around 0.1, but ACE-2 tails toward higher 


values consistent with urban and dust aerosol intrusion. 


BZ 


The TARFOX region has a noticeable mode at AOD around 0.2, 
but has significant spread of AOD values consistent with the 


varied aerosol constituents in that area. 


B. RECOMMENDATIONS 


As a result of this Seuay, the following 


recommendations are suggested: 


-Continue satellite retrieval validation with 
additional quality controlled data from these exercises, as 


well as future exercises. In situ data from ACE-1, TARFOX, 


and ACE-2 that have been well verified are increasingly 
becoming more available. ACE-Asia, which will be conducted 
off the west coast of Asia near Japan in 2000-2002 will 
provide additional regional information. This validation 
should include tests of individual components of the 
retrieval such as parameterization of the phase function, 


water vapor corrections, etc. 


- Assess the aerosol distribution model to see if a 
better fit can be obtained based on known aerosol 
distributions. Based on this and other recent studies the 
modes and standard deviations for the log-normal aerosol 
distribution can be tuned to specific regions or expected 


aerosol types. 


-Extend the S,, phase function parameterization 
technique to the next generation AVHRR with 1.6 um daytime 


channel. This wavelength is less sensitive to water vapor 


the the current channel 2. 


S13 


-Continue to improve on cloud screen analysis using 
more sophisticated techniques. Tests with static threshold 
values can be replaced with tests that determine thresholds 
dynamically from detected radiances. Thresholds can also be 


set by numerical model output such as for SST. 


-Improve sun glint removal process by incorporating 
real-time wind speed analysis. Wind speed as well as wave 
and swell direction can be obtained from other satellite 
sensors or from numerical models to get a more accurate sun 


glint approximation. 
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